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Abstract. We derive several identities for the Hurwitz and Rie- 
mann zeta functions, the Gamma function, and Dirichlet L- functions. 
They involve a sequence of polynomials otk(s) whose study was ini- 
tiated in [R]. The expansions given here are practical and can be 
used for the high precision evaluation of these functions, and for 
deriving formulas for special values. We also present a summation 
formula and use it to generalize a formula of Hasse. 



1. Introduction 
Let s G C and define aik(s) to be the Taylor coefficients given by: 

( - |at( ( '-" ) , "=i:»tf 1 w<i. (i.i) 

Then a (s) = 1, «i(s) — (s — l)/2, and, in jR], the following formula 
is proved: 



ak+i{s) 



1 A aj ( s )j(k + k 2 + s(2k + 2-j)) 



k(k + l)(fc + 2) (k-j + l)(k-j + 2) 



(1.2) 

The above recursion shows that, for k > 1, ak(s)/(s — 1) is a polynomial 
in s with positive rational coefficients, and allows one to obtain the 
bound [R]: 

(l + logCfc + l))!^ 1 
I"fc(s)|<c s , (1.3) 
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where 

c s = ^_M(| S | + 2)2N+\ (1.4) 

This bound was used in [R] in the derivation of the following formulas 
for the Riemann zeta and Gamma functions, valid for all s G C not a 
pole of the relevant function: 



oo 



rw = £^- (1-5) 







s + k 



T(s)^s + k 



For positive integer A: 



K ' k=0 j=l J 

where S(X,j) are Stirling numbers of the second kind. Finally, 



C(« + l) = =7-rj^a fc (s)* 1 ( S + A;), (1.? 



v y fc=0 

with + fc) the trigamma function. 

At the time of writing [Rj, the author also developed analogous for- 
mulas for Dirichlet L-functions via the Hurwitz zeta function, as well 
as a variety of additional formulas for the Riemann zeta function and 
the Gamma function, and we report on these here. We collect our main 
formulas in the theorem below. In equations (11. 9ft - (|1.16p . each stated 
formula is valid for any sGC that is not a pole of the lhs. 



Theorem 1.1. Let $lw > 0. Then, 

r(s + k + w) ' 



rw^T M rj #+fc ; . (i.9) 



fc=0 

Specializing to w — N + 1, a positive integer: 

T(s) = (N+lYN\Y- — , afcl f j -. (1.10) 

KJ V ; f^(s + k)(s + k + l)...(s + k + N) V ; 

K = 

For a > : 

oo 



cm = f^E^) r(g + fc + Q _ r - (i.n 
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For non-negative integer N : 

cw=£"-+^E ( . +t _ 1) .^ ifc -i +J o - (L12) 

For a > and A a non-negative integer: 

>/ i v-^ / n r(s + fc — 7 — i)r(a) . 

«« - *•«> = fw E £ rW^ cJA '* (L13) 



fc=0 jf=0 

where the coefficients c a (X,j) are defined by (15. 4p . 
H^e a/so /mve ; for positive integer A: 

EhC(.-A) = m E (, + t _ 2 )X*-A ' (1.14) 
a=i v ; k=o v 7 V ' 7 

where the b\ are given by (I6.12p . 

Let x be a non-trivial Dirichlet character for the modulus q. Then, 

L{s lX ) = a k {s)T{s + k-l)Y *(™) r (™/g) (1.15) 

V ,AJ q s T(s)^ KJ K ' ^ T (s + k + mq- 1 V ; 

^ w fc=0 m=l v ' ^ ' 

Furthermore, for A a non-negative integer: 

j oo g— 1 

L(s-A,x) = rf x A ^xMr(m/g) 

^ fc=0 m=l 

Er(s + - j - 1) 
o r( S + A; + m/g-j-l) Cm/9( 

(1.16) 

where c m / q (X,j) are defined by (15 .41) 

As a consequence of these formulas, we have several identities for 
special values. Two formulas for L(s,x) at non-positive integers s are 
summarized in the following theorem. 

Theorem 1.2. For non-negative integer r and a non-trivial Dirichlet 
character \ for the modulus q: 

L(~r, x) = r\q r Y f T I _ k )\ X(m)(m/q-l) . . . (m/q+k-r-1) 

k=0 ^ >' m=l 

(1.17) 
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and, for non-negative integer X: 

j=o \ J ' k=l KJ >' 

5-1 

x X(m)c m/q (\,j)(m/q- 1) ... (m/q- j + k). 

m=l 

(1.18) 

To compute a' fc (l) one can use either of the formulas from [Rj : 

k 

or ^ 

a ' k{i)= Yk\S, {x)kdx ' (L20) 

We also derive several other formulas for special values. For example, 

TV oo , x 

7 = £ I- W + 1) - wE g (1 . 21) 

m=l fc=l ' 

is shown, in Section HI to follow from fll.lOp . Equation (jl.2ip is known 
and first due to Kluyver |Klj . In [R], the author attributed his in- 
spiration to Kenter's short note on 7 |Kej . but the author has since 
discovered that Kluyver [Klj essentially had the same formulas as Ren- 
ter for 7, and also ( 11.21 j) . The above can also be viewed as a more 
precise form of the formula, often taken as the definition of 7, 

- 1 

7= lim > log(iV + l). (1.22) 

m=l 

Finally, we derive several interesting formulas that generalize a for- 
mula of Hasse. 

The alternating zeta function, or Dirichlet eta function, is defined 
by the Dirichlet series: 

17(a) :=E L -^ _ = C(a)(1_2l "' ) - (L23) 

n=l 

Theorem 1.3. For all s E C , and 9fts > — 1; 

F (*o + 1) to Vio (1 + exp(-x))^ dX J ^ (j + 1Y ■ 

(1.24) 
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Furthermore, for all s, s G C: 
r](s + s ) = 

E^-E^+mw^o+i-oeVtif 

m=0 ' 1=1 j=0 U ' 



where sim + 1, /) are Stirling numbers of the first kind. 



1.25) 



Notice that (jl.25p expresses rj(s), and hence ((s), in terms of its 
values at any collection of points so, sq — 1, sq — 2, 

The special case Sq = in the above theorem simplifies. Substituting 
t = exp(— x) into the integral of the first formula, 

f°° exp(-x(m + 1)) , , „. f 1 t m 

(m + 1) / — — -^rdx= m + 1 / — — -<ft 

1 J J 1 + exp -x m + 2 1 Vo l + t m+2 



t m+l 1 j 



we get 



Om+l 
Z 



1.26) 



00 1 m (_i+7 m ) 

m=0 j=0 VJ ' 

a formula first conjectured by Knopp and proved by Hasse [H] . later 
rediscovered by Sondow[S]. 

A similar formula holds for Dirichlet L-functions for a limited number 
of Dirichlet characters. For any non-trivial \ mod q < 5 we have: 

L( S + *o,x) = £^f £ < m + L 0£(*o + 1 - i, x) £ t-TT^' 

m=0 ' Z=0 j=0 V ' 

(1.28) 

valid for all s, So G C. The reason that we need to restrict here to q < 5 
has to do with the location of the g-th roots of unity relative to the 
point 1, as described in Section [8^1 



1.1. A summation formula. We give several other formulas related 
to the above. For example, we develop a formula that holds for all 
non-trivial Dirichlet characters for any modulus. To this end we use 
the following summation formula which we detail in Section [9j 

Let g, h be complex valued functions on the non-negative integers, 
say increasing at most exponentially, i.e. assume there exists c G R 
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such that g(n),h(n) = 0(exp(cn)). Define the m-th finite differences 
recursively: 

A°h(j) = h(j) 

A m h(j) = A m ~ 1 h(j + 1) — A m ' 1 h{j), m>l. (1.29) 



Then 



Let 



j2(-iyh( 3 + i)( m \ = (-i) m A m h(i). (i.3o) 

3=0 ^ ^ ' 



G(z):=J29(n)z n -\ (1.31) 

n=l 

We show that 

m 



Y,9{n)h{n)z n - 1 = J2G {m) (z)A m h(l)— } (1.32) 

n=l m=0 

holds in some disc centred on z = 0. 

In some applications, for example to the Hurwitz zeta function, one 
might prefer to write the n-th term as g(n)h(n)z n and start the sum at 
n — 0. Thus, for functions g, h on the non-negative integers, growing 
at most exponentially, 



Y,g{n)h{n)z n = Y J G {m \z)A m h(V)— v (1.33) 

n=0 m=0 

holds on some disc centred on z = 0, where 

oo 

G{z):=Y,9{n)z n . (1.34) 

n=0 

The analytic continuation of the rhs of (11.321) or fll.33f> to a given 
point z depends on the location of the singularities of G{w) or G(w) 
in relation to the point z, and also on the rate of growth of A m h(l) or 
A m h(0). In our applications, we substitute z = 1 after considering the 
analytic continuation of the function G or G. 

As a simple illustrative example, take g(n) = 1 if n = N and 
otherwise, so that G(z) = g(N)z N , and & m \z) = for all m > 
N. Substituting z = 1 we get, as a special case, Newton's forward 
difference formula: 

h(N) = /,(0) + iVAM0) + JV(JV 2 ,- 1} A> m + N(N ~ ^ N ~ 2) A 3 M0)+. ■ ■ 

Jl ' (1.35) 

As an application, we give the following identity. 
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Theorem 1.4. Let A > 0. If Us > 0, then 

Ms) = exp(-A) ^ A" ^ ( ~ ' _ (1.36) 

m=0 j=0 U 7 V Jy 

i/ie sum over m being uniformly convergent on compact subsets of Us > 
0. Furthermore, as m — )■ oo, 

y (-1)' (-l) m exp(-A) 

^ (j + l) s Ai(m-j)! A m (m+1 + A) s ' l ' J 

uniformly on compact subsets ofUs>0. 

Our summation formula is then applied to the Hurwitz zeta function 
in Section 19.11 and to Dirichlet L-functions in Section 19.21 

For the Hurwitz zeta function, we actually consider the alternating 
Hurwitz zeta defined by 

°° (_-\\n 
n=0 V ' 

where a > 0, and, in the above series, Us > 0. Our summation formula 
gives the following identities, as well as analytic continuation in the 
stated regions below. 

Theorem 1.5. For Us > —1 and s E C, 

C(s + s , a) = 



1 x 



x s ° exp(—x(m + a))(m + a + (a — 1) exp(— x)) 

ax 



r(«o + l)^ VVo (i + exp(-a;))™+ 1 



(1-39) 



For a// so, s 6 C, 

C(s + so, a) = J2 M~ E s ( m + !> OC*0»o + 1 - i, a) J] , ■ T \, • 
z — ' m! z — ' z — ' 7 + ar 

;i.4o) 



Let N + 1 fre a positive integer. Then, for all sq, s G C, 
C(s + s ,iV + l) = 



oo , s m m+l m (_l)j(m\ 



ml ' ' ' (7 + I 

m=0 Z=0 j=0 v 



;i.4i) 
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Notice that the second formula above is not a specialization of the 
first. The difference between the above two formulas is in the factor 
j '• + 1 vs j ' + a in the sum over j. 

For Dirichlet L-functions we have the following theorem, valid for 
any non-trivial Dirichlet character for any modulus q. 

Theorem 1.6. For3ts > —1, and s E C, 

L(s+s ,x)=y (m+1) / — ^ v . v dx > — ; 

(1.42) 

and, /or all s , s G C, 

, ^(-lr^ 1 , , 7W , A ^x(j + i)(7) 

L( a + ao ,x) = E L ^rE a ( m + 1 '0^o + i-0E fTTiV • 

m=0 ' 1=0 j=0 ^ ' 

(1.43) 

Notice that these expansions express L(s, x) i n terms of the eta func- 
tion 77(a) = C(s)(l-2^ s )! 

As a special case, set s = to get 

L ^ x) = £ ^tt E rTTTF ' (L44) 

m=0 j=0 VJ ' 

valid for any non-trivial Dirichlet character x, an d ah s G C. 

Formulas ( 11.240 - (11.280 and ( 11.390 - (11.440 converge uniformly for Sq, s 
in compact subsets of the regions in which the formulas are claimed 
to be true. As part of their derivation we give bounds for the finite 
differences given in each of the sums over j. 

Let t(x) be the Gauss sum 

9-1 

r(x) = Y,xUHj/q), (1-45) 

3=0 

where 

e(t) := exp(2vut). (1.46) 
We have the following estimates. 



Theorem 1.7. Let a > 0. Then, as m — > 00, 

(—iy I m\ log(m) s 
,_ U + a ) s \j) m a T(s J 

uniformly for s in compact subsets of C. In the event that s G Z and 
s < 0, then we interpret the ~ to mean equality when m > \s\. 
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Let x be a non-trivial Dirichlet character for the modulus q. Then 
as m/q 2 — > oo 



£* + u(7)~ 

e(-l/g)(l + e(-l/g)) m + X (-l)e(l/g)(l + e(l/g)) m ) 



j=0 



9 

< -icT\ (1.48) 



g l/2 <i 



and 

Furthermore, 



C q = |l + e(±l/g)| < 2. (1.49) 



j=0 v J / 



= lM 2 -+i e -^ x fooB(f(m + 2)), zfx(-l) = 1, 

q \-zsin(f(m + 2)), = -1. 

(1.50) 

as m/g 2 — >■ oo and m/q 4 0. In (I1.48P anc? (11.501) . we interpret ~ to 
mean equality whenever the rhs vanishes. 
We also have, for q fixed and as m — >• oo ; 

uniformly for s on compact subsets of 9fts > — M, where M is a non- 
negative integer, and with the implied constant in the <C also depending 
on q and M. 

These estimates say that the m-th finite differences of + a) s 
and of (— lYxU + 1)/ (J + l) s are exponentially smaller than the trivial 
bound obtained by using Y^™ (™) = 2 m . In the case of Dirichlet char- 
acters, the extra power of —1 is crucially needed to get to get such a 
small bound. 

We remark that Coffey, building on the author's paper [R] , has sep- 
arately arrived at formulas (11 .9p and (11.111) of Theorem 11.11 His appli- 
cation, however, is to the polygamma function and Stieltjes constants, 
whereas our focus is on formulas for L-functions and analytic aspects. 
Furthermore, his derivation is formal in that it does not contain any 
discussion on the convergence, or rate of convergence, of the sums in- 
volved in these formulas, and merely asserts analytic continuation. 
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Coffey uses, in his paper, the expansion 

'ln(x+l)V_ y> x k Bj z+k) 



x 



\x\ < 1, (1.52) 



k=0 



where the coefficients B^ +k ^ are generalized Bernoulli numbers given 



k 

by the generating function 



v 7 n=0 



(1.53) 



and remarks that 



(—-\) k n( s+k - 1 ) 

^M^.-D^L-i—. (1 .54) 

We regard the generating function (11. ip and the coefficients ak{s), to 
be more fundamental than the generating functions (1 1.53ft and (11.521) 
to the problem at hand. 

In this paper we focus our attention on deriving formulas for L- 
functions of degree 1. It seems, however, that some of our techniques 
can be applied to higher degree L-functions. We plan to explore higher 
degree L-functions, and to consider computational aspects of these for- 
mulas and variants, in a subsequent paper. 

2. Expansion for the Hurwitz zeta function 
The Hurwitz zeta function is defined, for 9fts > 1 and a > 0, by 

oo 

C(s,a) = J2(n + ar s . (2.1) 

n=0 

We have 

POO 00 

T(s)((s,a) = / X s ' 1 y^exp(— (n + a)x)dx 
Jo 



n=0 



-x^eM-ax)^ 



1 — exp(— x 

Substituting t — 1 — exp(— x), the above becomes, on applying (11.11) . 

C 1 dt 
T(s)C(s,a) = / (-logCl-f))- 1 ^-*)- 1 ^ 

Jo 1 

oo „i 

= Vor fc (s) / t s+k - 2 (l - tf-Ht. (2.3) 
k=o Jo 
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The rearranging of integration and summation over k is justified by 
bound f ll.3p which shows that, for given s, the sum over k converges 
uniformly for t in any closed subset of (0, 1). Now 

t s+k -\l ~ tr 1 = fi(s + k - 1, a) = l\ S l k ~ m f y (2-4) 

T(s + k + a — 1) 

Therefore, 

j-/ \ 1 , T(s + k - l)r(a) 

Note that (12.51) combined with the bound (II. 3ft for |«fc(s)| also provides 
the meromorphic continuation in s of ((s, a) to all of C, with the only 
pole being simple at s = 1, because the sum converges uniformly for s 
in compact sets away from its poles. For the latter fact, we also need 

r^*- 1 ' 1 (2.6) 



T(s + £; + a-l) k al 

as k — > oo which follows from Gauss' formula (I3.6p . also expressible as 

T(s + N + 1) ~ N S N\. (2.7) 

Apply the above with iV replaced by k, and s replaced, in the numerator 
of (12.61) . by s — 2, and, in the denominator, by s + a — 2. 

The poles of the numerator of the fc-th summand in (12.51) occur 
at s = —k + 1, —k, —k — 1, —k — 2, . . .. All of these, except for the 
pole at s — 1 when k = 0, are cancelled by the zeros of l/r(s) at 
s = 0, —1, —2, . . ., thus the only pole occurs at s — 1. 

If we let, above, a = 1 and use T(k + s — l)/r(s + k) = l/(s + k — l) 
we recover (jl.fip . And, if we take a = N + 1 to be a positive integer, 
we get 



rOO^(a + *-l)...(a + *-l + JV)- 



i.e. 



«') = E"- + ^)E (, +t - 1 ),°t' ) +t -i + iv) - (2 - 9) 

We can increase the rate of convergence of this sum by selecting iV to 
be larger. For example, in Table [1] we compare the precision achieved 
with s = 3, iV = 1, iV = 5, and iV = 20, iV = 100, and k < K, for 
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various values of K. More specifically, we truncate the sum at k < K 
and denote the remainder by : 

R(K,N,s):=-^- Y . ^ , — . (2.10) 

In our table, we depict the relative remainder terms 

R(K,N,s):=R{K,N,s)/((s), (2.11) 

for s = 3. 



K 


R(K, 1) 


R(K, 5) 


R(K, 20) 


R(K, 100) 





0.09876701304 


0.003735350802 


8.885612538e-05 


8.054816065e-07 


1 


0.06410420585 


0.001259436003 


1.056632936e-05 


2.148135338e-08 


2 


0.04503966189 


0.0005029064811 


1.595623563e-06 


7.505774529e-10 


4 


0.02594760777 


0.0001135959722 


6.161196116e-08 


1.697998088e-12 


8 


0.01207251038 


1.310916493e-05 


3.577149365e-10 


4.694895745e-17 


16 


0.004670251882 


7.985644458e-07 


2.203486644e-13 


1.938746948e-24 


32 


0.001580885488 


2.957659659e-08 


1.468535646e-17 


1.343384486e-35 


64 


0.0004908542321 


7.867242541e-10 


1.559347771e-22 


3.318995397e-51 


128 


0.0001445094596 


1.715803026e-ll 


4.578210941e-28 


2.669903582e-71 


256 


4.117178337e-05 


3.342699708e-13 


6.093017067e-34 


3.364327536e-95 


512 


1.148739118e-05 


6.112998825e-15 


5.183860931e-40 


7.923516397e-122 


1024 


3.160289401e-06 


1.078325101e-16 


3.464372739e-46 


3.670183638e-150 


2048 


8.606974887e-07 


1.861800819e-18 


2.037483910e-52 


1.820228675e-179 



Table 1. A comparison of the remainder terms 
R(K, N, 3) for = 1,5, 20, 100 and various values of 



K. 



Out of curiousity, we compare with the traditional way of computing 
((s) using Euler-MacLaurin summation: 

- ( S + ^ _1 ) [ B K {{t})r°- K dt. (2.12) 

where Bk,Bk(x) denotes the Bernoulli numbers and polynomials [R2J. 
and Bk := B&(0). Table [2] depicts, for s = 3 and the same values of 
A^ and K as in Table HJ the relative remainder term R2(K, N, s) : = 
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R 2 (K, N,s)/((s), where 

R 2 (K, N,s):=-( S + *~ ^ B K ({t})r s - K dt, 
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(2.13) 



K 










o 


-0 2478610589 


-0 002999137452 


-5 069543603e-05 


4 138739872e-07 


1 


0.1680926274 


0.0003284920385 


1.298774753e-06 


-2.079699113e-09 


2 


-0.03988421573 


-4.270910495e-06 


-1.080516795e-09 


6.931868307e-14 


4 


0.02944139866 


1.659288259e-07 


2.695929470e-12 


-6.931313924e-18 


8 


0.08490189016 


1.233350305e-09 


8.380153455e-17 


-3.464912761e-25 


16 


43.97791945 


3.248852315e-12 


4.036386548e-24 


-4.339822180e-38 


32 


2.000506814e+ll 


1.894001979e-13 


8.399102960e-35 


-6.233469067e-60 


64 


5.607820702e+39 


3.383888667e-07 


1.661297238e-47 


-6.266500089e-95 


128 


2.810377072e+115 


3.618530655e+24 


1.065733432e-54 


-1.163548513e-146 


256 


8.493688329e+304 


3.878203997e+124 


1.480620677e-31 


-1.227926982e-212 


512 


3.186217020e+760 


1.702666495e+401 


5.605030930e+91 


-1.276702408e-268 


1024 


2.130486030e+1825 


1.530602415e+1108 


2.919002699e+490 


-1.548683921e-227 


2048 


6.062592955e+4262 


7.835225059e+2829 


4.672401590e+1595 


-5.541412864e+162 



Table 2. A comparison of the remainder terms 
R 2 (K, N,3) for N = 1,5,20,100 and various values of 
K. 



Formulas ( 12. 9j) and ( 12.12ft are different in a number of ways. First, 
our sum over k converges and truncating the series after K terms 
gives a better approximation as K — > oo, while the sum in the Euler- 
MacLaurin formula provides a divergent asymptotic expansion for ((s). 
Second, our formulas give the meromorphic continuation of ((s) and of 
C(s, a) to all of C with a simple pole with residue 1 at s = 1. On the 
other hand, the Euler-MacLaurin formula only gives, for given K, the 
meromorphic continuation of ((s) up to > —K + 1. 

One can show (see 2.2.5 of [R2j ) . that the remainder term in the 
Euler-MacLaurin formula satisfies, for a := 5is > — K + 1, 

2((K) \s + K-l\ \T(s + K-l)\ 



\R 2 (K,N,s)\< 

1 2V ' ' n ~ N^ 1 a + K-1 |r(s)|(27riV)^ 
On the other hand, writing the fc-th term in (I2.10p as 

a k {s)T{s + k-l)T{N +1) 
T(s)T(s + k + N) 
we have, for given s and K, that 

\R(K,N,s)\ = 0(N-°- K ), 



(2.14) 



(2.15) 



(2.16) 
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with the implied constant in the big-0 depending on s and K. Here 
we have used the asymptotic formula 

r(s + k + N) n^- 1 1 ' ; 

as — 7- oo, which can be proven using Gauss' formula ( 13. 6 j) with s 
replaced by s + k — 1 or using Stirling's formula. Therefore, in the 
A^-aspect, formula ( 12. 8 p is roughly comparable to the Euler-MacLaurin 
formula, though with different dependence on K and s. The main 
difference between the two is in the extra (2ir)~ K in (j2.14jh and the 
afe(s), which grows exponentially in \s\, in each term of R(K, N, s). 

Next we examine the behaviour of R(K, N, s) in the K-aspect. The 
fc-th term of (|2TT0|) is 

a k (s)T(s + k-l)T(N + l) 

T(s)T(s + k + N) ' { ' 

By fll.3p . this equals, for fixed s and N, 

0(k~ N - 2+t ) } (2.19) 

for any e > 0, because ajt(s) = 0(k~ 1+t ) and T(s + k — l)/T(s + k+N) = 
0(k~ 1 ). Summing over k > K + 1 thus gives 

R{K, N, s) = 0(K~ N ~ 1+e ). (2.20) 

with the implied constant depending on s and N. 

3. A RELATED EXPANSION FOR THE GAMMA FUNCTION 

Let $lw, > and consider 

POO 

T(s)w- S = / x s ~ x e~ wx dx 







/ (-log(l-t)/*) a -V- 1 (l-t) w_1 d* 
Jo 

oo „i 
u-n JO 



t S+k -\l-t) W ~ l dt 



k=0 

oo 



Therefore, 



> Qifc(g) , r. (3.1, 

T{s + k + w) 



rM-THV^M. (3.2) 
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By (II. 3D . this series converges uniformly for all s in a given compact 
set away from its poles and therefore provides the meromorphic contin- 
uation of r(s). Specializing to w = N + 1, where N is a non-negative 
integer, we get 

r(s) = {n + iym V - — — — , a \\ } . — ; — — . (3.3) 

KJ v ; ^(s + k)(s + k + l)...(s + k + N) v ; 

fc=0 

Notice that, as iV — » oo, the k = term dominates. We can see this 
by writing the k-th term of the sum as: 

a k (s) 1 s(s + 1) . . . (s + k - 1) 



(s + k) s(s + 1) . . . (s + N) (s + N + 1) . . . (s + iV + jfe) ' 



(3.4) 



Now, s/(s + N + 1) = 0(1/N), with the implied constant depending 
on s. Furthermore, if N is sufficiently large in comparison to \dts\, then 
| (s + j)/(s + N + j + 1) | < 1 for all positive integers j. Thus the terms 
with k > 1 contribute 



^EfeM (3-5) 



(iV + l)*iV! lA|q t (s 
s(s + l)...(s + N)N f^\s + k 

to (13. 3j) . Bound ( 11.3j) shows that the sum over converges, hence the 
above is 0(1/N) times the k = term. We thus get Gauss' formula 
for T(s): 

iV s iV! 

r(s) = lim — -. (3.6) 

V ^ N^oc S (S + 1) ...(S + N) V ^ 

Here we have also used (N + l) s ~ A^ s as iV — > oo. Thus, the k = 
term in formula (I3.3P connects Euler's formula for the Gamma function, 
i.e. our starting point, to that of Gauss. 

4. Formula for Euler's constant 

In equation (13. 3p . subtract 1/s from both sides and then let s — > 0. 
On the lhs, we get —7, i.e. the negative of Euler's constant (here, we 
can define —7 to be the constant term in the Laurent series about s = 
of T(s)). On the rhs, the 1/s will cancel with a portion of the k = 
term. More precisely, expand (N + l) s as: 

(N + l) s = l + log(iV + l) s + CV). (4.1) 

Furthermore, the k = term, with the N\ but not the (iV + l) s , equals 

t^n, (4.2) 



s(s + 1) . . . (s + N) ' s + m \m 

\ 1 \ j m=0 
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the latter from the partial fraction expansion (I6.ip . Therefore, multi- 
plying f )4.ip . with (I3.3p . using the above for k = 0, collecting terms, 
and taking the limit as s — > 0, we have 



m=l x 7 fc=l 



k(k + l) ...(k + N)' 

(4.3) 

Now, by the binomial theorem, 

m=l v 7 

Substituting it = 1 + 1, expanding (m^ — 1)/(m — 1) = 1 + w + . . . u^ -1 , 
and integrating termwise gives 

E^rf!)= 1 + 1 /2 + ...l/iV. (4.5) 

Hence 

v 



m=l x 



7 . S I- M , + 1 )-,,g^)_. (, 6) 

5. HURWITZ ZETA FUNCTION SHIFTED 

Next we describe a formula for £(,§ — A, (3) where A is positive integer 
and a > 0. Begin with 



poo 00 

T(s)C(s - A,a) = / x 8 ' 1 + a) A exp(-(n + 

■™ n=0 



a)x)dx 



„1 oo 

/ (- log(l - t)) 3 - 1 YVn + a) A (l - t) n+a ^dt. 

Jo n=0 



(5.1) 



We can express 



^(n + a^l-t)™^- 1 (5.2) 

n=0 

as a rational function in t by starting with the A = case, Yl^Loi^ ~ 
t) n+a ~ l = (1— t) a_1 /t, and repeatedly multiplying by 1— t and applying 
—d/dt. We can thus prove, inductively, that 

f> + a)\l - t)^ = (1 - t)- 1 £ (5.3) 

n=0 j=0 
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where c a (0, 0) = 1, c a (0, 1) = 0, and 

c a (A + = (a - j - l)c a {X,j)+jc a (X,j - 1), (5.4) 

for A > 1, < j < A + 1 (for the recursion, we also set c a (A, A + 1) = 0. 
We need not worry about defining c a (A, j — 1) when j = because of 
the factor of j that appears in front of the c a (A, j — 1)). As usual, we 
write (-log(l -t))*- 1 = ET=o a k( s ) tk+8 ~\ and § et 

C(* - A, a) = w ^ afc (.) ^ _ _ c a (A, j). (5.5) 

When a = 1 we have 

Cl (A,j) = (-l) A+J j!5(A,j) (5.6) 
and (15. 5p can be seen as a hybrid of f ll.7p and (12. 5p . 

6. A LINEAR COMBINATION OF ( THAT CONVERGES QUICKLY 

Let A be a positive integer. In this section we derive a formula that 
expresses a linear combination of ((s — 1), ... , ((s — A) in a series of a 
nature similar to (I2.9p . 

We will exploit the fact that the same coefficients ak{s) appear 
in (II .7p independent of A, and that the denominators in (I1.7P are par- 
ticularly simple. By taking linear combinations of r(s)£(s — A) we can 
develop a sum whose terms converge more rapidly. We do so by using 
the partial fraction expansion: 

1 ^ + ^ + ...^^, (6.1) 



z(z — 1) ... (z — m) z z — 1 z — m 

where 

-l) m - 1 (-l) m ~ l fm 



ai = W^iy.=^^{i)- (6 - 2) 

The formula for a\ can be derived by considering the residue of both 
sides of (16.11) at z = I. 

Therefore, letting A be a positive integer, we wish to find a linear 
combination: 

A=l V ; k=0 A=l j=l J 



such that 

A 



g(-l)^ S (A, J b^ = glQ: i 1 ), (6.4) 
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A 


j 


c a (A, j) 








1 


1 





a — 1 


1 


1 


1 


2 





(a — l) 2 


2 


1 


2a - 3 


2 


2 


2 


3 





(a-l) 3 


3 


1 


3a 2 - 9a + 7 


3 


2 


6a - 12 


3 


3 


6 


4 





(a - l) 4 


4 


1 


4a 3 - 18a 2 + 28a - 15 


4 


2 


12a 2 - 48a + 50 


4 


3 


24a - 60 


4 


4 


24 


5 





(a - l) 5 


5 


1 


5 a 4 _ 3Q a 3 _|_ 7Q a 2 _ 75 a + 3i 


5 


2 


20a 3 - 120a 2 + 250a - 180 


5 


3 


60a 2 - 300a + 390 


5 


4 


120a - 360 


5 


5 


120 


6 


o 




6 


1 


6a 5 - 45a 4 + 140a 3 - 225a 2 + 186a - 63 


6 


2 


30a 4 - 240a 3 + 750a 2 - 1080a + 602 


6 


3 


120a 3 - 900a 2 + 2340a - 2100 


6 


4 


360a 2 - 2160a + 3360 


6 


5 


720a - 2520 


6 


6 


720 


7 





(a - I)"' 


7 


1 


7a 6 - 63a 5 + 245a 4 - 525a 3 + 651a 2 - 441a + 127 


7 


2 


42a 5 - 420a 4 + 1750a 3 - 3780a 2 + 4214a - 1932 


7 


3 


210a 4 - 2100a 3 + 8190a 2 - 14700a + 10206 


7 


4 


840a 3 - 7560a 2 + 23520a - 25200 


7 


5 


2520a 2 - 17640a + 31920 


7 


6 


5040a - 20160 


7 


7 


5040 



Table 3. A table of c a (X,j). 



for each 1 < j < A. Doing so would then result, by (16. ip . in the 
formula: 

tw-»^E ( , +t . 2) ,t +t -A-D - (-) 
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We now determine b\ so that (16.31) holds. Notice that (16.41) . for 
1 < j < A, can be written as the matrix equation 



DS 



/6i\ 



WJ 



:-i) 



A-l 



(A-l)! 



( { 7) \ 

-(V) 

V(-i) s:i (i:!)/ 



where S is the upper triangular matrix 
and D is the diagonal matrix 



D 



2! 



3! 



\ 



A!/ 



Thus, 



b 2 



W/ 



■1) 



A-l 



-l n-l 



(A-r 



with 



/ 1! 



LT 1 



1/2! 



-(V) 



1/3! 



1/Mj 



and 



5- 1 = ((-i) i+ M^)) 



AxA ' 



(6.6) 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



(6.11) 



where s(j,i) are the Stirling numbers of the first kind. Therefore, the 
solution to equation (I6.4p is given by: 



h _ (~1) A+A A sjj. A) (A - 1 



(6.12) 



thus yielding (16. 5p . 
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7. Expansion for Dirichlet L-functions 

Let x be a non-trivial Dirichlet character for the modulus q. Let 
L(s, x) be the Dirichlet L-function 

oo 

L(s, X ):=J2x(n)n- s . (7.1) 
n=l 

The defining sum is absolutely convergent in 3ts > 1 and, if x is n °t 
the trivial character mod q, conditionally convergent for 9fts > 0. Re- 
arranging the sum, for 3?s > 1, we can write it as 

1 9-1 

— ^2x(m)((s,m/q). (7.2) 

m=l 

Applying (12. 5ft and changing order of summation gives 

L(s,x) = > afc(s)r(s + fc-l) > — , -. 7.3 

r(s) <f ^ ^ r(s + fc + m/q - 1) 

Similarly, for any positive integer q, 

r W cw = 1 £ a, W r (s + 1 - 1) £ r( , + ^g } . (7.4) 

* fc=0 m=l v ' ^ 7 

Notice that the sum over m in (17.41) runs up to m < g rather than g — 1 
since, without the factor xi m )i we cannot ignore the m = q term. 

Next, if we substitute s — 1 into (17.31) then only the k — term 
contributes because ck^(1) = if k > 1, and r(s + — l)/r(s + k + 
m/q — 1) does not have a pole at s — 1 when k > 1. 

To deal with the = term, consider, for a > 0, the Laurent series 
about s = 1 

r(s-l) 1 7 + ^(a) 



T(s + a-l) r(a)(s-l) T(a) 

where 



+ 0(s-l), (7.5) 



:= ^V- (7.6) 
r(a) 

Substituting this, with a = m/q, into the = term of ( 17. 3ft . summing 
over m and using x( m ) = 0, so that the terms involving l/(s — 1) 
and 7 sum to zero, and then letting s — 1, we get the well known (see, 
for example, Proposition 10.2.5 of [C]) formula 

L(1 , x) = zlg xW ^M. (7 . 7) 
9 t^x r(m/g) 
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Next, let r be a non- negative integer. Consider the residue of (I7.3p at 
s = —r which involves, on the rhs, the terms with k — 0, . . . , r + 1. 
Using the fact that the residue at s = —r of T(s + k — 1) is equal to 
(-l) r - fe+1 /(r -k + 1)!, we have that 



L( _ r r) _ rW y M-r) \- X (m)T(m/q) 

L(r ' X) " r!g ^ (r + l-ifc)! ^Tim/q + k-r-iy {T8) 



We can simplify slightly using Y{m/q)/Y{m/q — r + k — 1) = (m/q — 
1) . . . (m/q — r + k — 1) for k < r + 1 and 1 for k = r + 1. The latter 
case does not contribute when summed against x( m )> an d we get 



L(~r, x) = r\q r ^2 ( r + i_ 53 X(m)(m/q-l) . . . (m/q+k-r-l) 

k=0 '' m=l 

(7.9) 

For example, when r = the above gives 

q-l q-l 

£(0, x) = - 53 X(m)(m/q - 1) = — 53 X(m)m. (7.10) 

m=l m=l 



The coefficients «&(— r) that are needed for a given r in (17. 9p are easily 
computed using, for instance, the recurrence (11.21) . 

We can write a similar formula for £(— r), though, without the x{ m ), 
we cannot ignore the k = r + 1 term, and, furthermore, the sum over 
m needs to be up to m < q: 



C (-r) = r \q r £ ( n^4^T " ^ • • " ^ + fc " r ~ ^ 

fc=0 ^ m=l 

(7.11) 

The innermost summand above is taken to equal 1 if k — r + 1. For 
example, if q = 1, we get 



C(-r) = (-l) r r!a r+1 (-r). (7.12) 



In, Section 7 of [R], the author showed that a r+ i(— r) = B r+1 / (r + 1)!. 
Thus, the above gives Euler's formula £(— r) = (— l) r £> r+ i/(r + 1). 
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Next, we develop a second formula for L(s, x) evaluated at negative 
integers. Using ( 15. 5 p we have, for positive integer A: 

-y j oo q-l 

L(s-X,x) = T(s) s —\ Yl ak ^ Yl xijn)Y{m/q) 

^ ' V k=0 m=l 



ET(s + k — j — It 
Q 1 (s + k + m/q — j — 1) 



(7.13) 



Substituting s = 1, we get contributions from the terms k = 0, and 
from j = k with j = 1, . . . , A. The latter is accounted for by the pole 
of T(s + k — j — 1) balancing out against the zero of ak(s) at s = 1. 
Thus, 

j=o V J ' fc=i u y ' 
x ^ x ( m ) Cm /<z( A >-7')( m /9- X ) • • • ( m /^ - J + 

m=l 

(7.14) 



8. General Dirichlet series 

Let 

b(n) 



n 

71=1 



be a Dirichlet series, assumed to be absolutely convergent for 3fts > 01. 

Substitute formula (13.11) . with w = n into the Dirichlet series, and 
change order of summation over n and k: 



-. oo oo 



a fc (s)r(s + A;)(n-l)! 



oo oo 



1 / n , / .T(s + k)(n — 1)\ 

- r^E-WEW r(. + t + n) • (8 - 2) 

v ; fc=0 n=l v ' 

Changing order of the sums can be justified for !Rs > a±, but we do 
not do so here, because our emphasis below will mainly be on the first 
equality, rather than the second. However, before turning our attention 
to the first equality, we note that, if b{n) is identically equal to 1 for 
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n > N + 1 and otherwise, then the inner sum over n in the last line 
above equals 

T(s + k-l)N\ 



+ k + N) 

This can be seen by taking the identity 

T(s + fc)r(n) 



/\ OO oo 

n=N+l n=N+l 



T(s + k + n) 



.3) 



U) 



and summing the geometric series on the lhs. Thus, in this case, equa- 
tion (18.21) reduces to (12. 9p . The fact that we can simplify the sum over 
n is what makes (18. 2p particularly useful in the case of the Riemann 
and Hurwitz zeta functions, and Dirichlet L-functions. 

Next we explore the first equality of (18. 21) . Using the partial fraction 
expansion (16. ip . we have 

(n-1)! (-I)" (n-l\ (g5) 

T(s + k + n) ^j + Jfe + ml m )' 

Therefore, 

^^EwE^wE^kf; 1 )- (8 - 6) 

v ' n=l k=0 m=0 1 1 \ / 

But bound (11.31) allows us to rearrange the two inner sums over k and 
m: 

E Ms) E ( n ~ 1 )= E(-ir f n_1> )E - J ¥ £ ^- 

(8.7) 

The resulting inner sum over k equals 

-1 poo 

;-log(l -^/ty-H'+^dt = / a; s - 1 exp(-s)(l -exp(-x)) m dz 

- ™g£&G> (8 - 8) 

which we see using (11.11) . then changing variable t = 1 — exp(— x), 
expanding (1 — exp(— x)) m using the binomial theorem, and integrating 
termwise. We thus have the identity 

1 V ak ^ = V I m \ (R Q\ 

T{s)^s + k + m ■ [ -> 
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Notice that, while in the integrals in ( 18. 8p we require 9fts > —m, both 
sides of the above equation are entire functions of s, and hence, by 
analytic continuation, equality holds for all s. 
Hence, 

oo n-1 , 1 x m (_-l\j(m\ 

*M=E*oE(-ir("; 'Era(r- («<» 

n =l m =o v / j=0 u y 

Since the inner double sum is just another way to express 1/n 8 , for- 
mula ( 18.1 Op holds for all s for which the defining Dirichlet series con- 
verges conditionally. 

We can also arrive at the above more directly by starting with the 
binomial coefficient identity 



(8.11) 



for j > 0, and n > 1. To prove this identity, consider the coefficient of 
z n ~ l in (z - l)"" 1 YZ =j O zm = (-l) n "V(l - z) n -i- 2 . 
Thus 

E^EW^OtT)^. <-> 

because only the j = n — 1 terms survives with an extra (—1)" intro- 
duced from the inner sum. Substituting into the Dirichlet series for 
L(s), and rearranging the inner double sum gives (18.101) . 

8.1. Heuristic manipulations. Next we proceed heuristically, but 
we will justify our resulting formulas in certain cases in the next section. 

Changing order of summation, and ignoring for the moment the fact 
that the resulting sums over n diverge unless b(n) is rapidly decreasing, 
gives 

n_1 (-i) m 00 m (-l) j ( m ) 

L ^ = J2-^H - 1) •••(---) E ttti^- (8 ' 13) 

m=0 ' n=l j=0 ^ ' 

We interpret the typically divergent sum over n as being equal to the 
m-th derivative, evaluated at z — 1, of the analytic continuation, if it 
exists, of the series 

oo 

f{z):=Y,Kn)z n -\ (8.14) 

71=1 
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thus suggesting 

m=0 '' ' j=0 VJ ' 

For example, let L(s) = ((s)(l - 2 /2 s ) = 1 - 2~ s + 3" 3 - A~ s + . . . 
be the alternating zeta function. We consider here the alternating zeta 
function, rather than ( itself, because the pole of ((s) = 1 creates a 
slight complication that we wish to circumvent. The factor (1 — 2/2 s ) 
has a zero at s = 1, balancing the pole of ((s). Thus, C( s )(l ~ 2/2 s ) 
extends to an entire function. Then f(z) = 1 — z + z 2 — z 3 + . . . = 
1/(1 + z), and / {m) (l) = (-l) m m!/2 m+1 , hence (Km reduces to 

C(.)(l-2/2*) = £ — (8.16) 

m=0 j=0 U ' 

i.e. to Hasse's formula described in the introduction. 
We can also, in (18.131) . expand 

^ . m+1 

(n— 1) . . . (n—m) = —n(n—l) . . . (n—m) = — s(m+l, Z)r?/. (8.17) 
n n 1=0 

Here, we have multiplied and divided, in the first equality, by n so as 
to have, in the second equality, a polynomial in n rather than in n — 1. 
Substitute the rhs of (I8.17P into (I8.13p . and rearrange the sum over n 
and / to get 

m=0 ' Z=0 n=l j=0 VJ ' 

We interpret the typically divergent sum Y^=i K n ) nl1 as being equal 
to L(l — I), suggesting the following interpolation formula for L(s): 

°° (_-\\m m + 1 m (-1)3 ( m ) 

L (») = E i 7^ L E s ( m + - £ 7-TTF' (8 ' 19) 

m=0 ' Z=0 j=0 VJ ' 

Now, s(m + 1, 0) = for all m > 0, and thus we can start the sum at 
I = 1 unless L(s) has a pole at s — 1. 

Consider again the case L(s) = C( s )(l — 2/2 s ). In this case, for / > 1, 

L(l - /) = C(l - 0(1 - 2') = (-1)^(1 - 2 l )B t /l, (8.20) 

and the sum over m in (I8.19P equals 

m+1 

^(-l)'- 1 ^ - 2 l )s{m + 1, l)Bi/l. (8.21) 
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Comparing with (18. 16ft suggests that the above equals 

(-l) m m!/2 m+1 . (8.22) 

In the next section we will rigorously prove the above manipulations in 
a few cases, including the alternating zeta function, and thus yielding 
equality between (lOTj) and f[Q2j) . 

Note that we can also consider the above manipulations for the 
Dirichlet series with coefficients b(n)/n s ° for some so G C. Writing 

L(s + s ) = jr b Ml. (8 .23) 

n=l 

and letting 

n=l 

we anticipate (in certain cases) two formulas: 

+ ( s, 5 ) 

m=0 ' j=0 VJ ' 

where L^ m '(so,l) is the m-th derivative with respect to z of L(sq,z) 
evaluated at z = 1, and 

00 (_-\\m m + l m ( — \)3 ( m ) 

l(s+s ) = E [ -Jr E < m +^ iW'o+i-i) E TTTT^- (8 ' 26) 

m=0 ' 2=0 i=0 ^ ' 

8.2. Justifying the manipulations. We are able to prove the formu- 
las of the previous section, specifically equations (I8.25P and f !8.26|) for 
the alternating zeta and alternating Hurwitz zeta functions, and Dirich- 
let L-functions L(s, x) where is non-trivial character for the modulus 
q < 5. In Section l§T2l we develop a related summation formula that can 
be used to give interesting formulas valid for all q. 

For higher degree L-functions, for example, degree 2 L-functions 
associated to cusp forms of given weight, level, and character, it appears 
from numerical experiments that formula f!8.26|) holds but with some 
adjustment- additional convergence producing terms, arising from a 
smoothed approximate functional equation, seem to be needed. We 
will revisit the issue of higher degree L-functions in a future paper. 

We start by modifying equation (I8.10p by introducing an extra pa- 
rameter z. Let 

b(n) 



L( S ,,):=E^ n_1 - ( 8 - 27 ) 



n=l 
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Assume that there exists a real number r such that 

6(n) < n r , (8.28) 

so that the coefficients b(n) do not grow too quickly, and the series 
defining L(s,z) converges for all s G C and \z\ < 1. Then, as before, 
replacing l/n s by ( 18. 12ft . we rigorously have 



- l\ ™ ( i)'(7) 



CXJ II— 1 / \ 

n — 1 <m — V / 



x , ^ ti + 1)" ' 

n=l m=0 x x j=0 VJ y 



(8.29) 



Now, for \z\ < 1/2, the sum over n and m converges absolutely. That is 
because ( I1.47P shows that the sum over j is bounded, in the m-aspect, 
by an amount which is generously 0(1), with the implied constant 
uniform for s in compact subsets of C. Therefore, the absolute value 
of the rhs of the above is, on using Ylm=o ("m) = 2 n_1 , 

oo 

^^Tl&Hl^r 1 . (8.30) 

71=1 

This converges for all \z\ < 1/2 because b(n) is assumed to grow at 
most polynomially in n. Therefore, for \z\ < 1/2, and s 6 C, we have 
on rearranging the sums over m and n and pulling out a z m /m\ from 
the sum over n: 

°° ( — y\ m °° m (— l) j ( m ) 

m=0 n=m+l j=0 ^ J ' 

(8.31) 

Our next goal is to obtain the analytic continuation in z of the sum 
over n and to substitute z — 1. In general, this cannot be done, but in 
some examples we are able to do so. 

8.3. Alternating zeta function. The alternating zeta function is de- 
fined to be 

rt«):=£M— ■ (8-32) 

n=l 

We wish to allow shifts, so we let b(n) = (— l) n_1 /n so for some Sq G C. 
Thus we are considering 

00 (—\)n-l ~7i-l 

v (s + s ,z):=J2-^7 ( 8 - 33 ) 

71=1 

This can also be regarded as a special case of the Lerch zeta function. 
We have introduced the parameter sq since we wish to examine the 
general expansions of the form considered in (18.261) . 
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In this case, the sum over n in (18.311) reduces to 

{ -^—(n-l)...(n-m)z n - m - 1 . (8.34) 

n=m+l 

Now 

T = J- r X "" leXp( ~ X) ix (8 35) 

^ r(s )y l + 2exp(-x) ' 

which we can prove by expanding the denominator as a geometric series 
and integrating termwise. The above integral converges and is analytic 
in both variables for IRso > and z ^ (— 00,— 1]. Differentiating the 
lhs termwise m times within its disc of convergence gives (18.341) . Hence 
differentiating the rhs m times yields 

V ^ '- (n - 1) . . . (n - m)z n - rn - 1 

n=m+l 

'-l) m m\ f°° x 30 - 1 exp(-x(m + 1)) ^ 



r(s ) 7 (l + zexp(-s)) m + 1 

Also note that we can bound these derivatives using knowledge of the 
radius of convergence of the Taylor series, about the point w = z of 
the function t](sq,w). Consider the Taylor expansion 



OG 



' ml 

m=0 

Here 77 ^ refers to the m-th derivative of 77 with respect to the second 
variable. Assume that 9ftz > —1, so that the nearest singularity of the 
rhs of (I8.35P is the point —1. Thus, the above converges absolutely for 
all \w — z\ < |1 + z\, in particular, at w — if \z\ < |1 + z\. The latter 
holds when §lz > — 1/2 (as can be seen by drawing z and 1 + z). 

Thus, returning to flOTj) (with L(s) = r)(s)), for Rz > -1/2, Rs > 
0, s e C, 

- r (so ) 2w* U (i + ,ex P (-,))^ dx )^( j + iy 

(8.38) 

(recall, for the purpose of convergence, that the sum over j is O s (l)). 
Therefore, substituting z = 1 we get 

„(, + so) = j- f ( r ^ eM - x{m + 1)} d X ) f 

W + °°) T(s ) ^ Uo (1 + «p(-x))^ ^ (j + 1)- ' 

(8.39) 
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thus justifying ( 18. 25ft in the case of L(s) = T](s). For any given so 
satisfying 9fts > 0, the rhs above converges uniformly in s on com- 
pact subsets of C, thus also giving the analytic continuation of the 
alternating zeta function. 

We can extend the validity of (18.391) in the s aspect by integrating 
by parts 

00 X s0 ' 1 exp(-x(m + 1)) , (m + 1) f°° x s ° exp(-x(m + 1)) , 

dx = / — -, rr — -rrdx. 



(1 + exp(-x)) m+1 s Jo (1 + exp(-a;)) m+2 

(8.40) 

While the integral on the lhs is convergent for 3ftso > 0, the integral on 
the right converges for 5Js > — 1. Therefore, for 9fts > —1 and s G C, 

v(s + s ) = 

f (m + 1) ( r xS0 exp( " x(m + 1)} **) T 

T ^ + !) ho Vo (1 + exp(-x))^ dX ) (j + 1)' ■ 

(8.41) 

For example, substituting s = gives Hasse's formula (I1.27P as de- 
scribed in the introduction. 

Because the m-th term decreases exponentially in m for given sq 
(with 3fts > — 1), and uniformly for s in compacts subsets of C, equa- 
tion ( I8.4ip also gives the analytic continuation of rj. 

The above also shows that, for any s G C and Jfeo > — 1, the limit, as 
z — > 1 of the rhs of (I8.33P exists and is equal to rj(s + s ) as expressed, 
for example, in ( I8.4ip . Therefore, returning to ( I8.34p and expanding 
using ( I8.17P we get 

i](s + s ) = 

°° I nm^ 1 m (-l)j( m ) 

m=0 ' 1=1 j=0 U ' 

(8.42) 

While we initially assumed Jfeo > — 1, we can extend ( I8.41|) to 3ftso > 
—M, where M is any positive integer, by repeatedly integrating by 
parts. Therefore, equation (18.421) holds for all sq, s G C. 

8.4. Dirichlet L-functions. Let \ be a non-trivial Dirichlet character 
for the modulus q. Define 

%^) = E^ Z "" 1 . 1^1 < 1- (8-43) 

n=l 
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Now, for Kso > 0, 

which we can see by expanding the denominator as a geometric series 
and using the periodicity mod q of x- Note that, for given x, the 
singularity of the integrand at z = exp(x) is removable because the 
numerator vanishes at z — exp(rc). Therefore, the rhs is analytic in 
some neighbourhood, depending on q, of z = 1. 

As in the discussion concerning 77, we can get an estimate for the 
size of the m-th derivative with respect to z of L(s , z, x), evaluated at 
z = 1 by considering the singularities of the rhs of (18.441) which occur 
along rays emanating outward from the non-one q-th roots of unity, i.e. 
at z — rexp(2irij/q), r > 1, 1 < j < q. 

We require that the distance from the point 1 to the nearest singular- 
ity, exp(±27ri/g) be greater than 1 so that the sum analogous to (18.381) 
continues analytically to z — 1. This occurs when q < 5 (at q = 6, 
exp(27ri/6) has a distance of 1 from the point 1, and for q > 6 the 
closest root is even closer). Thus, for any non-trivial x mod q < 5 we 
have: 

L(s + s , x) = £ [ -Jr £ s(m + 1, l)L(s + 1 - i) £ 

m=0 ' 1=0 j=0 \J ' 

(8.45) 

valid for all s, s G C. 



9. A SUMMATION FORMULA 

Next we consider a generalization of (I8.12p . The same binomial 
coefficient identity which gave that formula for l/n s can be used for 
any function h(n) on the positive integers: 

h (n) = J2(-iyhu + 1) £(-ir ( n ~ m l ) (7) , (9.1) 

j=0 m=j ^ \ J / 

because only the j — n — 1 term survives. Therefore rearranging sums, 

n— 1 / _. \ m 

m=0 ^ ' j'=0 
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Let g, h : N — > C be two functions on the positive integers. Apply 
the above identity for h(n) to the power series 

oo 

Y,g{n)h{n)z n -\ (9.3) 

n=l 

and then rearrange the sum over m and n, pulling out al/m! and z m 
from the sum over m: 

71=1 

= E^T E 9(n)(n-l)...(n-m)z n - m - 1 

m=0 n=m+l 

x£(-w + i) 7). 

3=0 ^ J ' 

(9.4) 

The latter step of rearranging the sums over m and n can be justified 
for sufficiently small z assuming modest bounds on g, h. For example, 
assuming that g(n) and h(n) are 0(n r ) for some real number r, and 
using trivial bounds on the sums of binomial coefficients, the sum over 
m and n converges absolutely for \z\ < 1/4. By trivial, we mean using 
Y™ =0 (J) = 2 m , and J2Z=o = ^ ■ 0nce can even assume that 
g, h grow exponentially, i.e. 0(exp(cn)) for some csR, and still get 
an identity valid on some disc centred on z = 0. 

We can view the sum over n as being the m-th derivative with respect 
to z of the power series 

oo 

G(z):=Y,9(n)z n -\ (9.5) 

n=l 

Thus equation (19.41) can be written as: 

oo oo / N m m / x 

S^)Mn)^ 1 = E G(m) ^ i ^E(- 1 ) i/l 0" + l)(7)- ( 9 - 6 ) 

n=l m=0 ' j=0 

Note that, in the previous sections, with the choice /i(n) = l/n s , and 
the asymptotic formula (OTjl . so that E^o(- 1 ) i (7)/(j + !) s < s ^ 
we were thus allowed to initially take \z\ < 1/2. The second ingredient 
for i] and L(s,x), with ^f(n) = (— l) n ~ l /n s ° or x(n)/n s °, was to obtain 
the analytic continuation of the rhs of (19. 4p to a neighbourhood of the 
point z = 1. We achieved this by expressing the appropriate G(z) in 
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closed form in ( 18. 35ft and ( I8.44p from which we were able to deduce the 
rate of decay of G^ m \z). 
The inner sum over j, 

j2(-i) j Kj+i) r y (9.7) 

can be interpreted in terms of the finite difference operator. Define 

aMj) = h(j + i)-h(j) 

A m h(j) = A m ~ 1 h(j + 1) - A m ~ l h(j). (9.8) 



Then <$1§ is equal to 
hence 



-l) m A m /i(l), (9.9) 



^)p(n)A(n)« n - 1 = ^ G^(z)A m /i(l)— . (9.10) 

n=l m=0 

Note that, because our application is to Dirichlet series, we prefer, 
above, to write the n-the term of the lhs as g{n)h(n)z n ~ 1 . However, in 
other applications, as to the Hurwitz zeta function below, one might 
prefer to write the n-th term as g(n)h(n)z n and start the sum at n — 0. 
Thus, as a variant, for functions g, h on the non-negative integers 



oo r 

Z 



V g(n)h(n)z n = V G im) (z)A m h(0) — , (9.11) 
^— ' m! 

n=0 m=0 



where 



G(z) := X>(n)s». (9.12) 

n=0 

The range of validity in z of (19.101) and (19. lip depends on the rates of 
decay of g(n), and A m h(0) as discussed above, for example being valid 
on some disc centred at if g(n) and h(n) grow at most exponentially 
in n. The analytic continuation of the rhs to a given point z depends 
further on the location of the singularities of G(w) or of G(w) in relation 
to the point z. 

To illustrate our summation formula, let A > 0, and take g(n) = 
(-1)"-^", h(n) = n~ s \- n . Then G(z) = X/(1+Xz), and G( m )(l)(-l) m /m! 
A m+1 /(1 + A m+1 ). Formally, this choice of g and h gives Amore's gen- 
eralization of Hasse's identity [A"] : 



+ A ) m+1 ^ + X ) SAJ+ 



(9.13) 
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In order to justify the above, one also needs to obtain a rate of decay 
for the sum over j (depending on A), and this is done in [X] . 

Many variants are possible. For example, let A > 0, g(n) = (— l) n_1 
X n /(n - 1)!, h(n) = n" s A" n (n - 1)!. Then G(z) = exp(-Xz), and 
G <ym \l) = (— A) m exp(— A). Therefore, substituting into our summation 
formula, and simplifying the factorials and binomial coefficient that 
appear, we get 

oo m f-lV 

7]{s) = exp(-A) Y \ m y - -f— ? -. (9.14) 

m=0 j=0 KJ ' v J ' 

We justify this formula for 9fts > 0. Note that, while h(n) grows very 
quickly with n, to justify the above formula what we need, for the 
purposes of convergence and rearranging sums as in (\9A\i , is an estimate 
for the size of the sum 

y 015) 

fe(j + iyXf(m-j)V { ■ ] 

For > 0, this sum is equal to 

Note that, as m — > oo, the main contribution to the above integral 
arises from smaller x, which we can see by breaking the integral into, 
say, f Q + Jj 00 , and on the second integral, breaking the sum into < 
j < m/2 and m/2 < j < m. We thus truncate the integral at x — 1, 
and then reverse the sum over j to get, as m — > oo, 

j2 



X m T(s) 



j (j + iy\j( m -j)\ 

[\^eMHrn+l)x)y { - X)j T {jX) - ( 9 ' 17 ) 

The resulting sum in the integrand is the truncation of the Maclaurin 
series for exp(y), with y = — Aexp(x), i.e. 

^ (-^expW =exp( _ Aexp(x)) _ £ (-XYeMjx) (9lg) 

j=0 ] ' j=m+l 3 ' 

An analysis shows that, for given A > and as m — > oo, that the 
latter sum does not contribute to the main asymptotics of the inte- 
gral. Therefore, dropping that sum, and then extending the domain of 
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integration to oo, we have 

-IV' 



i \m roo 

' _s-l 



(-l) m exp(-A) 



x s exp(— A exp(x) — (m + l)x)dx 

(9.19) 



A m (m + 1 + A) 5 ' 

with all our estimates above uniform on compact subsets of 9fts > 0, 
thus justifying equation (19.141) in that half-plane. The last ~ can be 
derived by substituting t = (m + l)x into the integral and expanding, 
for t = 0(m + 1), 

A exp(i/(m + 1)) = A + \t/(m + 1) + 0(\t 2 /m 2 ). (9.20) 

A more careful analysis shows for a given A, that to get a good 
approximation (i.e. within a specified small relative error e), for a 
given A, we need m ^> A, with the implied constant depending on s 
and e. 

Formula (ll.36p becomes interesting as we take A large, because the 
terms in that formula are, for m sufficiently large and including the 
exp(— A) that appears in front of the sum, 

(9.21) 

{m+l + X) s K J 

Thus, while (— l) m /(m + 1 + A) s decays at roughly the same rate as 
terms in the Dirichlet series defining rj(s), our estimate shows that the 
terms in ( 19. 14ft are eventually, for m ^> s A, exponentially small. Thus, 
for A large, there is a transition zone, as we sum over m, in which we 
quickly achieve roughly 2A/log(10) digits accuracy for r)(s) before the 
sum continues its leisurely convergence to i](s). 



9.1. Application to the alternating Hurwitz zeta function. We 

define, for a, 3?s > 0, the alternating Hurwitz zeta function to be 

(*(s,a) :=Y / \ . (9.22) 
v ' ^ (n + a) s v ' 

n=0 v ' 



IDENTITIES HURWITZ AND L-FUNCTIONS 35 

Setting, in (ETTTj) . g(n) = (-l) n /(n + a) s ° and h{n) = \/{n + a) s , we 
have 



■l) n 



n + a) s ° 

— — / X s "- 1 V(-2) n exp(-(n + a)a;)da 

dx, (9.23) 



n=0 

-oo 

n=0 

1 f°° x 30 ' 1 exp(-ax 



r(s ) 7 l + zexp(-x) 



and so 



^ (n + a) s ° (n + a) s 

n=0 v ' v ' 

1 f^ m ( r x-Q- 1 exp(-x(m + a)) i \ A (~1) J (7) 
r («o)^ Uo (l + ^exp(-x))™+ 1 x )j^(j + a y 

(9.24) 

We remark that, here, the function is equal to $(—z, sq, a), where 
$ is the Lerch zeta function. 

A similar analysis of ( I9.23|) to that of ( I8.35P can be carried out and 
shows that fl9~24j) holds for s 6 C, 5?s > and z <£ (-oo, -1]. The 
asymptotic formula (jl.47p provides a uniform bound for the sum over 
J- 

Thus, substituting z — 1 gives 

(*(s + s , a) = 

1 " / /- ^- 1 exp(-x(m + q)) \ A (-1) J (?) 
r(«o) ^ U (1 + «p(-s))«+i V ^ (J + a)* • 

(9.25) 

We can extend the range of validity to 3fts > ~~ M by repeatedly inte- 
grating by parts M times. For example, for 3fts > — 1 and s G C: 



C*(s + s , a) 



x s ° exp(— x(m + a))(m + a + (a — 1) exp(— x)) 



r( So + l)^ Wo (l + exp(-x)) m + 1 

: ')-(;;•) 

(j + a) s ' 



3=0 
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Substituting s = 0, the above integral can easily be evaluated to equal 
l/2 m+1 , thus specializing to another formula of Hasse. 
We also have, as in Section IHI3| for all s , s G C, 



OO / -, \m "H-l 



m f_l\jfm\ 

C(s + s , a) = J2 E < m + !> OC*0»o + i-u)E 7 ■ \ [i ■ 

m! U U (j + a) 

(9.27) 

We give yet another formula for the alternating Hurwitz zeta func- 
tion, when a = N + 1 is positive integer. Taking, in our summation 
formula fl£3D, 9( n ) = if n < iV, and (-l) n - N ~ 1 /n s ° if n > iV + 1, 
and h{n) = l/n s , we have: 

00 00 ' -zY 1 ^ 1 



G(z) = Y,9{n)z n - X = {-I) 



N \ ^ 



1 JV+1 

= r., ^(-(*+i)*) dx . (9,8) 

r(s ) Jo 1 + ^exp(-x) 

This gives the analytic continuation of Cr(z) to the point z = 1 where 
G has radius of convergence 2. Therefore, substituting z — 1 and 
proceeding as before 

C*(s + so,iV + l) = 

E ^r- E •(»• + 1. DCi" + 1 - » + 1) E 77TTF' 

m=0 ' Z=0 i=0 VJ ' 

(9.29) 

9.2. Dirichlet L-functions again. Next we give formulas that ex- 
press L(s,x), for any non-trivial character \ for the modulus q (i.e. 
without restriction on q), in terms of the Riemann zeta function, or, 
more precisely, in terms of the alternating zeta function. 

We let g(n) = (— l) n_1 /n s ° and h{n) = (—l) n ~ l x{n)/n a . Therefore, 
G(z) here is identical to the function in ( I8.35P and the whole discussion 
of (I8.32p regarding the function G(z) carries through. We also require 
a bound for the sum over j, which we will show beats the trivial bound 
exponentially. 

Let 

C q := |l + e(l/g)|. (9.30) 
Assume that 5is > — M where M is a non-negative integer. We will 
prove, for q fixed and as m — > oo, 

±^±^( m )«m M C- (9.31) 
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uniformly for s on compact subsets of 9fts > — M (with the implied 
constant in the <C also depending on q and M). But C q < 2, hence, 
dividing by 2 m shows that Yl'jLo '{j+iy (7) decreases exponentially 
fast as m — > oo. It is crucial for this bound, and also for getting a nice 
rate of decay on G^ rn \z)/m\, that we include here the factor (— l) n_1 
in g(n) and h(n). 
Thus: 

30 exp(-x(m + 1)) \ ^ xU + 1) (7) 



^ Vio (1 + exp(-i-))™+ 2 



(9.32) 



and 



^ ^(-ir^ 1 , 1 .s / , A fx(j+i)0 

L( S + , , X ) = £ £ a (m + 1, /)i7(ao + !-*)£ (? + 1 y " 

m=0 ' Z=0 j=0 ^ ' 

(9.33) 

Equation H9.32I) converges for Jfeo > —1, and s G C, while (I9.33f) 
converges for all so, s 6 C, uniformly for so, s on compact subsets of C, 
thus also giving the analytic continuation of L(s, x)- 



9.3. A bound on Y^™ xU + 1) (7) / U + l)* - I n ^ s sec ^ on we use the 
notation 

e(t) := exp{2mt). (9.34) 
We have the following identity 

X(j + 1) /m N 



E 



„ (j + i)'Vj 

<? — 1 m -i 1 \ / \ 1 

'x(a)e((a-l)l/q), (9.35) 



lyy e(-jl/q) ( m \y^ 
H 1=0 7 =0 u y v J 7 a=l 



<7 

which is easily verified by noting that, for given a, that only the terms 
j with j + 1 = a mod q survive the sum over I. Each a, above, thus 
contributes 

j + l — a mod m 

Summing over all 1 < a < g thus gives the lhs of (19.351) . 
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We first examine the s = case because we can obtain very precise 
information in that case. Identity (I9.35P reads 

m , x 

£*(, + u (7 

3=0 VJ/ 

= i£(l + eH/g)r^x(o)e((a-l)Z/g), (9.37) 

^ 1=0 a=l 

Notice, crucially, the sum over a on the rhs vanishes when 1 = 0. It is 
for this reason that we are able to beat the trivial bound. 

Furthermore, |l + e(— l/q) \ < 2 for 1 < I < q — 1, with the largest val- 
ues being when I = 1, q— 1. Also note that, above, the sum over a can be 
expressed in terms of the Gauss sum. When I = 1 the sum over a equals 
e(— l/q)r(x), and when / = q — 1 the sum equals e(l/g)x(— 1)t(x). 

Therefore, for fixed q, 

^ ( e (-l/g)(l + e(-l/g)) m + X (-l)e(l/g)(l + e(l/g)) m ) 

where C g = |1 + e(±l/g)| < 2. Note that we have also used |r(x)| < 
g 1//2 , with equality when \ is primitive. Dividing by 2 m , the rhs decays 
exponentially in m. Therefore, because G(z) here is given by f!8.35j) and 
has Taylor series with radius of convergence 2 at z = 1, we can sub- 
stitute z = 1 in our summation formula thus yielding equations (19.321) 
and f)9.33p . when s = 0. 

We can be a bit more precise about the rate of decay of the middle 
expression of (I9.38p . Expand 

1 + e (t) = 2(1 + nit - n 2 t 2 + ^vr 3 t 3 + 0(t 4 )). (9.39) 
The bracketed term can be written as 



exp(mt - n 2 t 2 /2 + 0(t 4 )). 



(9.40) 
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Therefore, 

^ (e(-l/g)(l + e(-l/q)r + *(-l)e(l/ ff )(l + e(l/g)) m ) 

= ^) 2V ¥( 1+0 ( 1 /' 2 )»(e-i( m+2 »+x(-l)ei (m+2) ) 

= l(x) 2 -+ 1 e ~^ (1+0(1/<?2)) x J cos (^ m + 2 ))' if = 

g ^ \-isin(J(m + 2)), ifx(-l) = -l. 

(9.41) 

Note that in the case that x is a rea l primitive character, the above 
simplifies further to 

2™+ 1 c -^ (1+0(1/g2)) x (cos(^(m + 2)), if = 1, 

g 1 ^ 6 9 |sin(|(m + 2)), if x (-l) = -1. 



(9.42) 

We can similarly bound the contribution from the terms 2 < I < q — 2 
in (I9.37p . For simplicity assume that % is primitive so that the inner 
sum in (I9.37|) can be written in terms of t(x) and is, in absolute value, 
equal to g 1 / 2 or 0, depending whether gcd(/, q) = 1 or not. Furthermore, 
amongst 2 < / < q — 2, the maximum |1 + e(l/q)\ m occurs when I = 2 
or q — 2. Thus the contribution from 2 < I < q — 2 is bounded by 

2V / V 2 ^ (1+0(1/92)) , (9.43) 

i.e. the rate of decay in m for the remainder term is roughly 4 times 
larger. 

We have thus shown, for primitive \i that 



m / \ 

!>■+<;) 



= 2m+ i e -^(i+o(i/ g 2 )) x | cos(^(m + 2)), if x(-l) = 1, 

q ^ \-zsin(^(m + 2)), ifx(-l) = -l. 

+i?( X ,m), (9.44) 
where the remainder R(x, m) satisfies 

i?( X ,m)<2V / V 2 ^ (1+0(1/92)) . (9.45) 

With respect to uniform asymptotics allowing both q and m to grow, 
equation (I9.44p provides an asymptotic formula for the sum on the left 
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as m/q 2 — > oo and m/q 4 — )■ 0, 




cos(^(m + 2)), 
-zsin(£(m + 2)) 



if = 1, 

if = 



-1, 



(9.46) 



while the second line of (19.38P provides an asymptotic formula valid 
in the larger region m/q 2 — > oo. Note that we do no need to assume 
that x is primitive in the above asymptotic formula- that assumption 
was made to simplify the discussion regarding R(x,m), but a similar 
analysis applies. 

We interpret the above asymptotic formula to be an equality in the 
event that the rhs vanishes. There are two cases to consider. If 1) = 
— 1, then the rhs vanishes iff m — —2 mod q. On the other hand, the 
terms j and m — j on the lhs cancel each other out, assuming that 
j 7^ m — j , because the binomial coefficients match, and, furthermore, 
whenm = -2 mod q, x(m-j + l) = x(-2-j + l) = x(-l)x(j + l) = 
~x{j + !)■ Special care is needed when m is even and j = m/2 because 
then there is just one binomial coefficient, and not a pair. But, if 
j = m/2 then x(j + 1) — x(( m + 2)/2). However, m + 2 is assumed, 
here, to be a multiple of q, hence gcd((m + 2)/2, q) > 1 (because q > 2 
when x is non-trivial), and so x(( m + 2)/2) = 0. Thus, xU + 1) — 0, 
and the middle binomial coefficient does not contribute to the sum. 

The case x{~ 1) — 1 is a bit more complicated, but the idea is the 
same, namely to show that the terms j and m — j cancel. We have 
cos(^(m + 2)) = iff m + 2 = rq/2 where r is an odd integer, i.e. 
2m + 4 = rq. Thus g is even because r is odd, and, in fact, divisible 
by at least 4 because there are no even primitive characters for the 
modulus 2. Hence m = —2 mod q/2. Because r is odd, we can write 
it as r = 2r + 1. Now, x{ m ~ J + 1) — x((^ r o + 1)9/2 — j — 1) = 
x(g/2 — j — 1) = x(j + 1 — 9/2), the last step because x(~ 1) = 1- 

Next, factor g into distinct prime powers, q = with 
pi = 2, and let x(n) = Now, for 1 < / < fc, 

"1 ^/e 

+ 1 - 9/2) = XJ>iU + 1)) because g/2 = mod pf J for I > 1. 

Thus, to prove here that x( m ~ J ' + 1) = — x0 + 1)> i- e - that the terms 
j and m — j cancel, we need to show that 



Xaft0* + 1 -?/2) 



X 2 ft (j + 1) 



(9.47) 
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for all j and any primitive character mod 

Note that we can assume that j + 1 is odd because q/2 is even. Thus, 
if j + 1 is even then both Xtfi (j + 1 — q/2) and Xzh (j + 1) are equal 
to 0. 

Write q = 2^(2g + l), i-e. write the odd part of q as 2q + 1. Thus, 

x 2 h (3 + 1 - g/2) = x 2 ft (j + 1 - 2*" 1 ). 

Now, if /?i = 2, then ^(j + 1) = — XaU — 1) and we are done. Next, 
if (3\ > 3 then any character for the modulus 2 /?1 > 8 can be written, 
for odd n, as 

^W = e(^- + ^) (9.48) 

where < < 2, < // < 2/ 31 " 2 , and 

n = (-l)^' mod 2^, (9.49) 

with < ^ < 2, < z/ < 2 /3l ~ 2 See Chapters 4-5 of [D]. Furthermore, 
X2-31 is primitive implies that // is odd. 
We claim that 

n - 2 ft-1 = n^ 1 - 1 + 1) mod 2 pl (9.50) 

which holds because n is assumed to be odd, and 

5 2^- 3 = 2 ft-i + i mod 2^, (9.51) 

which can be proven inductively by repeatedly squaring. 
Therefore 

Xaft (n - 2*- 1 ) = xHx^- 1 + 1) = xHx(5 2 " 1_3 ). (9.52) 

But 

(5 2 " 1 " 3 ) = e(^|^) = exp(vr^') = -1, (9.53) 

the last step because fjf is odd when x 2 ft is primitive. 

We have therefore shown that if: x is primitive, x{~~ 1) = 1> an d 
cos(^(m + 2))= 0, then x(j + 1) — — x( m ~~ j + 1) f° r an J- Finally, we 
need to address, as before, the case where m is even and j = m/2, i.e. 
when the middle binomial coefficient is not part of pair. But, in that 
case, x(j + 1) = x(( m + 2)/2) = x(rg/2) which equals because g is 
divisible by 4 and hence q/2 is even. 

9.3.1. Bounding the sum in general. Assume for now that 9fts > 0. 
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We consider the sum over j in (19.351) . which we can write as 

m 



e{-jl/q) ( m 

, (1 u ■ s V./ 



r( s ) 



x s_1 (l + exp(-x)e(-//g)) m exp(-x)cfe. (9.54) 



Similar to the s = case, one can show that the I = I, q — 1 terms 
dominate (the / = term in (I9.35j) vanishes), and we have, for fixed q 
and as m — > oo, 



^ (i + i) s U 



3=0 



r(x) 



e(—l/q) / X s (1 + exp(-x)e(-l/q)) m exp(-x)dx 



qT(s 

+ x(-l)e(l/g) / ar' -1 (l + exp(-x)e(l/g)) m exp(-x)rfx ). 



(9.55) 

Taking absolute value and using |l+exp(— x)e(— l/q)\ < C q gives (jl.5ip 
Repeated integration by parts of (I9.54p . say M times, allows us to 
obtain the same sort of bound for > —M, at a cost of, on integrating 
by parts, introducing extra powers of m into the integrand, thus giving 
bound (Oil) . 

This rate of decay, beating 2 m by an exponentially small amount, 
justifies formulas f l9.32|) and (I9.33p . 

10. An estimate for E£o£^(7) 
For a > 0, we derive an asymptotic formula for 



as Tfi — y oo. The above equals, for m + 3?s > 0, 

s s_1 exp(— ax)(l — exp(—x)) m dx = 

(- \og(v/m)) s ~ 1 v a - 1 (l - v/m) m dv. (10.2) 



1 f°° 

r( a ) 
i 



m a T(s) 

Here we have substituted v = mexp(-x). Now, we can write 

logim/vy- 1 = logm s_1 ((l - logw/logm)^ 1 -1 + 1), (10.3) 
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so that ( 110. 2ft equals 

Mm)- 1 v a-i {l _ v/m)mdv + 

((1 - log v/ log m) 3 - 1 - Ijv^il- v/m) m dv. 



m a T(s) J 



m a T(s) J 

(10.4) 

It is a simple exercise to show that the first integral above is asymptot- 
ically equal to T(a) as m — > oo by breaking up the integral into, say, 

1/4 

rm ' , rm , 

Jo +J m i/ 4 , and using: 

(1 — v/m) m = exp(mlog(l — v/m)) 

= exp(-,-|--|--...), \v\<m 

= exp(-u)(l + 0(l/m 1/2 )), \v\<m 1/4 . (10.5) 

We can also show that the second integral in (110.41) does not contribute 
to the main term. First, undo our change of variable and express the 
integral as 



oo 



((x/\ogm) s 1 — l)(mexp(— x)) a (l — exp(— x)) m dx. 



o 



We need to show that this tends to as m — > oo, and uniformly for s in 
compact subsets of C. To do so, break the above into three integrals: 

»log(m)— log(m) 1 / 2 







oo 



log (m) +log (m) 1 ' 2 

+ / . (10.6) 

log(m)— log(m) 1 / 2 J log(m)+log(m) 1 / 2 

In the middle integral, we have [xj log m) s ~ l — 1 = 0(log(m) -1 / 2 ), with 
the implied constant uniform for s in compact subsets of C. Hence, 
pulling this estimate out, and using, as before, 

(mexp(-x)) a (l - exp(-x)) m cfe ~ T(a) (10.7) 



o 

we have 

"log(m)+log(m) 1 / 2 



((x/\ogm) s 1 - l)(mexp(-x)) a (l - exp(-x)) m dx 

log(m)— log(m) 1 / 2 

0(log(m)- 1/2 ). (10.8) 
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Next we bound the first integral, using 

(1 - exp(-:r)) m = 

exp(— m(exp(— x) + exp(— 2x)/2 + exp(— 3rr)/3 + •••)) 

< exp(— mexp(— x)). (10.9) 

Note, further, that (1 — exp(— x)) m increases monotonically for x > 0, 
which, combined with the above estimate gives 

(1 - exp(-:r)) m < exp(- exp(log(m) 1/2 )), (10.10) 

when < x < log(m) — log(m) 1 / 2 . This shows that (1 — exp(— x)) m 
decreases much faster than any power of m, in particular m a . Further- 
more, if 5fe > 1, then (rr/logm)*" 1 — 1 is bounded over the interval 
[0,log(m) - log(m) 1/2 ]. If Iks < 1, then let 3?s > —M, where M is 
a non-negative integer. For m > M, we can steal M + 1 powers of 
(1 — exp(— x)) = x + 0(x 2 ), for |x| < 1, from (1 — exp(— x)) m , with 
little effect on that factor, so that (x/ logm) 5 " 1 — 1)(1 — exp(— x)) M+1 
is, on the interval [0, 1] of size 0(log(m) _CT+1 + 1), where o = 3fts. Fur- 
thermore, the latter O estimate also holds for [xj logm) 41 " 1 — 1) once 
x > 1, if < 1. The integrand of the second integral is thus handily 
dominated by the last factor of size 0(exp(— exp(log(m) 1//2 ))) and the 
second integral tends to 0, uniformly for s in compact subsets of C. 

Finally, to bound the third integral, consider the factor [xj logm) s_1 — 
1. If 3Rs < 1 then this is bounded for x > log(m) + log(m) 1 / 2 . And 
if 5Rs > 1 then this factor is 0(x a ~ 1 / log(m) CT_1 ). Using the trivial 
bound (1 — exp(— x)) m < 1, we thus have an upper bound for the third 
integral, when < 1, of 

a f°° , exp(-alog(m) 1 / 2 ) 
< m a / exp{-ax)dx = — ^ — - (10.11) 

Jlog(m)+log(m) 1 / 2 a 

and, when 9fts > 1, of 

< i — r=i / x*- 1 cxp(-ax)dx. (10.12) 

log(m) CT Jlog^+logim) 1 / 2 

Substituting u = x — log(m) — log(m) 1 / 2 , then pulling out (log(m) + 
log(m) 1 / 2 ) " -1 from the integral, and finally using 1 + y < e y , with 
y = u/{\og(m) + log(m) 1//2 ), the above is bounded by 

. 1/2 (log(m) + log(m) 1 / 2 ) CT " 1 
exp(-alog(m) 1/2 ) v toV '- — - toV / — 

r < < (^-!) 

x / exp(w(-a + - — — - , 1/2 ))du 
Jo (log(m) + log{my/ 2 

<exp(-alog(m) 1/2 ), (10.13) 
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with the implied constant, for given a, uniform for s in compact subsets 
of C. 

We have thus shown, for given a, that the second integral in (110. 4p 
tends uniformly, for s in compact subsets of C, to 0, hence proving the 
asymptotic formula (jl.47p . 

Note that when seZ and s < 0, then the first line in ( 110. 2p . valid 
for m > —s, vanishes because 1/T(s) is then equal to 0. Therefore, 
the asymptotic formula (11.471) should be interpreted as an equality for 
such s and m. 
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